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Abstract. We present first results of a partial wave analysis of the diffractive reaction n~ Pb — > n~n + n~ Pb 
based on data from the COMPASS experiment taken during a pilot run in 2004 using a 190 GeV/c n~ beam on 
a lead target. The analysis was performed in the region of squared four-momentum transfer t' between 0. 1 and 
1.0 (GeV/c) 2 . The n~n+n~ final state shows a rich spectrum of well-known resonances. In addition a spin-exotic 
J PC = V + state with significant intensity was observed at 1.66 GeV/c 2 in the p(770)/r decay channel in natural 
parity exchange. The resonant nature of this state is manifest in the mass dependence of its phase difference to 
jpc _ j++ an( j 2-+ waves. The measured resonance parameters are consistent with the disputed ^i(1600). An 
outlook on the analyses of the much larger data set taken during 2008 and 2009 is given. 



1 Introduction 

The naive Constituent Quark Model (CQM) describes light 
mesons as bound color-singlet states of quarks and anti- 
quarks with flavors u, d, and s grouped into SU(3)fl av0 r 
multiple ts. In the CQM the spin J, parity P, and charge 
conjugation C of a meson are given by 

(1) J = \L-S\,...,L+S; P = (-1) L+1 andC = (-l) L+s 

where L is the relative orbital angular momentum of quark 
and antiquark and S = 0, 1 the total intrinsic spin of the 
qq' pair. In addition a meson is characterized by its isospin 
/ and G-parity which is defined as 

(2) G = (-1) I+L+S 

Both quantum numbers are conserved in strong interac- 
tions. 

Despite of its simplicity, the CQM works surprisingly 
well and explains a good part of the observed light me- 
son spectrum as well as some of the meson properties, al- 
though it does not make any assumptions about the nature 
of the binding force. Quantum ChromoDynamics (QCD) 
describes the strong interaction between colored quarks 
by the exchange of gluons, which are colored themselves. 
Even though it is still not clear, how the confinement of 
quarks and gluons into hadrons is generated by the QCD 
Lagrangian, its structure suggests that there are color-sing- 
let states beyond \qq). In particular one expects gluonic 
degrees of freedom, that reflect the non-Abelian charac- 
ter of QCD, to manifest themselves in the meson spec- 
trum. Mesonic states beyond the CQM are classified into 
hybrids, glueballs, and multi-quark states [1]. Hybrids are 
\q~qg) resonances with constituent glue that contributes to 
the quantum numbers of the hadron. The qq' pair is in a 
color-octet configuration which is neutralized in color by 
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the excited gluon. Glueballs are pure gluonic bound states 
\gg). Multi-quark states include tetraquarks and mesonic 
molecules. Most of these states will be hardly distinguish- 
able from ordinary \qq') states with the same J PC and will 
mix with them so that physically observable states are in 
general a linear combination of \qq') and additional basis 
states that is experimentally very difficult to disentangle. 
An unambiguous evidence for the existence of mesonic 
states beyond the CQM, as allowed by QCD, would be the 
discovery of exotic states with quantum numbers forbidden 
in the simple quark model. In the light-meson sector so- 
called spin-exotic states with quantum numbers J PC = , 
+ ~, 1~ + , 2 + ~, 3~ + , 4 + ~, ... are a particularly promising 
field of research. 

In lattice QCD [2] simulations the lightest glueball is 
predicted to have ordinary scalar quantum numbers J pc = 
++ and a mass of about 1.7 GeV/c 2 . A possible experi- 
mental glueball candidate is the / (1500) seen by the Crys- 
tal Barrel [3] and WA102 [4] experiments. The interpreta- 
tion of the data is, however, complicated by the mixing of 
the / (1500) with other states. 

The lowest mass hybrid, on the other hand, is predic- 
ted [5] to have exotic quantum numbers J PC = 1~ + so that 
it will not mix with ordinary resonances. It is expected to 
have a mass between 1.3 and 2.2 GeV/c 2 and in the flux- 
tube model it preferentially decays into Z7i(1235)tt and 
/i (1285) 7r [6]. Three experimental candidates for isovec- 
tor spin-exotic J pc = 1~ + states were found so far: The 
7Ti(1400) was seen in n~N -> n~r]N by the E852 [7] and 
VES [8] experiments. Crystal Barrel observed a 1~ + rj7T 
state in pn — > tTt^t] [9] and pp — > npTpri [10] Dalitz 
plot analyses. Another 1~ + state, the 7Ti(1600), was seen at 
higher mass by the E852 and VES experiments in 
pn [11,12,13], rfn [14,15], f x n [16,17], and b x n [18,17] 
decay modes in peripheral n~p interactions. It was also re- 
ported in pp — > binn [19]. The resonant nature of both 
1~ + states is still discussed controversially in the commu- 
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Fig. 1. Production of resonance X in diffractive scattering of the 
beam particle /?beam off the target particle h tmgst and its subsequent 
decay into n hadrons h\,...,h n . The interaction is mediated by 
exchange of a Reggeon R and leaves target particle intact. At 
high energies the Pomeron P is the dominant Regge-trajectory. 

nity [8,17]. In particular the observation of the 7Ti(1600) in 
the pn: decay channel is heavily disputed. A different anal- 
ysis of a larger E852 data set showed that an extension of 
the wave set, used to model the data, makes the ;ri(1600) 
signal disappear [20]. A third 1~ + state, the 7Ti(2000), was 
seen only by the E852 experiment in the fin [16] and 
bin [18] decay channel. 

The COMPASS experiment sets out to settle these and 
other issues in light-meson spectroscopy by studying me- 
son production in diffractive and central production reac- 
tions using high-energetic hadron beams on various fixed 
targets. 

Diffractive reactions are known to exhibit a rich spec- 
trum of produced states. In diffractive meson production 
the target particle remains intact and interacts with the in- 
cident beam particle via /-channel Reggeon exchange. In 
this process the beam particle /ibeam is excited to some res- 
onance X which then dissociates into the final state as il- 
lustrated in Fig. 1: 

(3) h he . dm + /itarget -» X + h[ mget and X -> hi + . . . + h„ 

The process can be characterized by two kinematic vari- 
ables: the square of the total center-of-mass energy, s = 
(Pbeam + ^target) 2 , an d the squared four-momentum transfer 
from the incoming beam particle to the outgoing system X, 
t = (Pbeam - Px) 2 ■ It is custom practice to use the variable 
t ' — \t\-\ t | m ; n instead of t, where 1 1 | min is the minimum value 
of \t\ allowed by kinematics for a given X invariant mass 
m X - The value of |f| min is small but larger than zero, be- 
cause of the additional longitudinal four-momentum trans- 
fer which is required by m x > mb eam . At high beam ener- 
gies the diffraction is dominated by Pomeron exchange [21] 
so that isospin and G-parity of the intermediate state X are 
the same as that of the beam particle. In addition the final 
state particles are produced mostly at small angles with re- 
spect to the beam direction which requires a high angular 
resolution and coverage of the detector. 

In central production reactions the mesonic state is pro- 
duced in Reggeon-Reggeon fusion, where both beam and 
target particle remain unaltered. This is illustrated in Fig. 2. 
These reactions will enhance scalar mesons and in par- 
ticular to the / (1500). It is also believed that in central 




Fig. 2. Central production of meson X in Reggeon-Reggeon fu- 
sion and its subsequent decay into n hadrons h\,...,h n . Both 
beam and target particles stay intact. 

production gluonic degrees of freedom are enriched which 
makes it a promising reaction for glueball searches. The 
produced resonance X carries only a small fraction x F « 
of the maximum possible longitudinal momentum /?™ ax in 
the center-of-mass frame, where x F = Pl/p™** ~ ^-PlI V*- 
The scattered beam hadron has x F ~ 1 and the target x F ~ 
-1. In fixed target geometry this means that the scattered 
beam particle appears as the leading hadron hf ast in the 
event, while the target recoil proton is slow. The centrally 
produced system is separated from the leading hadron and 
the recoil particle by rapidity gaps. 

2 Experimental Setup 

Finding mesonic states that do not fit into the CQM in the 
light-quark sector is an experimentally challenging task. 
There is a high density of broad overlapping states so that 
the signal can only be extracted from interference effects. 
This requires large data sets and a complete phase space 
coverage of the experiment. The COmmon Muon and Pro- 
ton Apparatus for Structure and Spectroscopy (COM- 
PASS) [22] is a fixed-target experiment at the CERN Super 
Proton Synchrotron (SPS) and is well suited for this task: It 
is a two-stage spectrometer built to measure at high beam 
intensities. The detector covers a wide range of scattering 
angles and particle momenta and has a high angular res- 
olution. In the target region scintillating fibers and planes 
of silicon microstrip detectors are used for beam definition 
and vertexing. The setup has two large dipole magnets with 
1.0 and 5.5 Tm bending power which both are surrounded 
by staggered tracking detectors with increasing granularity 
and resolution towards the beam axis. The two spectrome- 
ter stages are equipped with hadronic and electromagnetic 
calorimeters. This enables COMPASS to reconstruct final 
states with charged as well as neutral particles like n°, rj, 
and rj'. In addition the first stage features a Ring Imaging 
CHerenkov (RICH) detector able to separate charged final 
state pions and kaons in the momentum range between 5 
and 50 GeV/c. 

The M2 beamline that feeds the COMPASS experi- 
ment is very versatile and can deliver high-intensity sec- 
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ondary hadron as well as polarized tertiary muon beams 
which are produced at a production target by the incoming 
400 GeV/c proton beam from the SPS. The hadron beams 
can have a momentum of up to 300 GeV/c and a maximum 
intensity of 4 ■ 10 7 sec -1 . The positive hadron beam con- 
sists of 71.5 % p, 25.5 % 7T + , and 3.0 % K + , the negative 
beam has 96.0 % iC and 3.5 % K . 

The suitability of the COMPASS setup for meson spec- 
troscopy was studied in a short pilot run in 2004, where 
a 190 GeV/c n~ beam was shot onto a fixed lead target. 
Based on the experience from this run the spectrometer 
was upgraded in order to address the challenges of the 
hadron spectroscopy physics program. A Recoil Proton De- 
tector (RPD) was installed around a 40 cm long liquid hy- 
drogen target. The RPD measures the time-of-flight of the 
recoil protons using two barrels of scintillator slats and 
provides information for the trigger decision. The read- 
out electronics of the electromagnetic calorimeters was up- 
graded and the central part of the second calorimeter was 
equipped with 800 radiation-hard Shashlik blocks. In ad- 
dition the tracking close to the beam axis and the vertex 
definition were improved by adding high-resolution Pixel- 
GEM detectors in the spectrometer and cryogenic silicon 
microstrip detectors in the target region, respectively. At 
the same time the material budget in the beam region was 
reduced. Also the particle identification capabilities were 
enhanced by utilizing the RICH detector in the first spec- 
trometer stage and two ChErenkov Differential counters 
with Achromatic Ring focus (CEDAR) upstream of the tar- 
get which are able to identify the incoming beam particles. 

With this enhanced setup COMPASS took diffractive 
and central production data in 2008 and 2009 using 
190 GeV/c negative and positive hadron beams on liquid 
hydrogen, nickel, tungsten, and lead targets. The goal of 
collecting about ten times the available world statistics for 
both production reactions has been achieved. 



3 Diffractive Production of 



a 4 cm central hole rejected beam particles that were not 
crossing the target material. Additional lead-scintillator ve- 
to counters placed downstream of the target suppressed 
events with outgoing particles that were emitted under large 
angles so that they would fall out of the acceptance. A sys- 
tem of three scintillation detectors in the spectrometer, so- 
called beam killers, vetoed non-interacting beam particles. 
In order to enrich diffractive events, a scintillator disk with 
5 cm diameter was placed after the target and was used 
as a multiplicity counter, selecting events with at least two 
forward charged particles. In addition at least one cluster 
with minimum energy deposit of 6 GeV was required in 
the hadronic calorimeter in the second spectrometer stage. 

The 2004 data sample contains about 87.7 • 10 6 events 
from the diffractive trigger which were subject to further 
offline event selection. Events were required to have ex- 
actly one primary vertex with one incoming beam and three 
outgoing charged tracks. The primary vertex position was 
constraint to the target region and the three outgoing tracks 
were required to have a charge sum of - 1 . Diffractive events 
were enriched by an exclusivity cut. However, the spec- 
trometer setup was not able to measure the energy of the 
incident beam particle event-by-event, only the direction 
of the beam particle and thus the scattering angle 6 was de- 
termined precisely by a silicon microstrip beam telescope. 
In the 2004 setup also the target recoil particle was not de- 
tected. For small momentum transfer, large beam energy, 
and assuming specific target and recoil masses the beam 
energy Ebeam an d the momentum transfer f can be calcu- 
lated in good approximation from the well-measured quan- 
tities nix, Ex, and 6. For the analysis the t' region between 
0.1 and 1.0 (GeV/c) 2 was selected, where the 7T!(1600) 
was reported in the past. In this f region the beam par- 
ticles were assumed to scatter off quasi-free nucleons in- 
side the Pb nuclei of the target. The validity of this as- 
sumption is supported by the roughly exponential falling 
f'-spectrum shown in Fig. 3 which has a slope parameter 
compatible with the nucleon radius. Exclusive n~7i + n~ pro- 
duction was selected by requiring the reconstructed beam 



During the 2004 pilot run COMPASS was able to record 
within a few days a data sample of the diffractive reaction 

(4) n~ + Pb -> n-jT + 7T- + Pb 

with competitive statistics using a 190 GeV/c n~ beam 
on a total 3 mm thick lead target. The n~7T + 7T~ final state 
was chosen for a first analysis [23,24,25], because the con- 
troversial spin-exotic 7T!(1600) has been observed in this 
channel. 



3.1 Data Sample and Event Selection 

The trigger for diffractive events selected one incoming 
and at least two outgoing charged particles. In 2004 it con- 
sisted of several components. Incoming beam particles were 
selected by a coincidence of two scintillator discs with 
5 cm diameter that were positioned upstream of the tar- 
get and centered about the beam axis. A veto system with 
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Fig. 3. Spectrum of the squared four-momentum transfer /' = |f| - 
|f| min in logarithmic scale. The data exhibit a roughly exponential 
behavior. 
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Fig. 4. Distribution of the beam energy £b ea m calculated from m x , 
E x , and 9 for the f range between 0.1 and 1.0 (GeV/c) 2 . The 
exclusivity cut is indicated by the vertical lines. 
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Fig. 5. n~n + n~ invariant mass distribution of the selected data 
sample for/' 6 [0.1, 1.0] (GeV/c) 2 . 



energy £beam t0 be maximum 4 GeV off the mean beam 
energy of 189 GeV (see Fig. 4). About 420000 events in 
the mass range between 0.5 and 2.5 GeV/c 2 pass all of the 
above selection cuts. 

In the 2004 pilot run there was no particle identifica- 
tion for the incoming beam particles. The non-pionic com- 
ponent of the beam of about 5 % consists mainly of K . 
Their decays and diffractive reactions constitute part of the 
background in the analysis. Also for the final state no par- 
ticle identification was applied. 



3.2 Partial-Wave Analysis 

Figure 5 shows the 7r~7r + 7r~ invariant mass distribution of 
the selected data sample. It exhibits clear structures in the 
mass region of the well-known resonances ai(1260), 
«2(1320), and 7T2(1670). In order to find and disentangle 
the various resonances in the data, a partial-wave analysis 
(PWA) was performed. 
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Fig. 6. A partial wave in the isobar model. The diffractively pro- 
duced X~ with quantum numbers J pc M e decays first into an iso- 
bar and a bachelor n~ . The isobar has spin 5 and a relative orbital 
angular momentum L with respect to the tfCacheior' The i s °bar then 
decays in a second step into n + iC . 



The PWA employed in this analysis is based on two 
basic assumptions. We assume that the total cross section 
factorizes into a resonance and a recoil vertex and we use 
the isobar model to describe the X~ decay. This by defi- 
nition neglects any final state interaction of the outgoing 
pions or the isobar with the target as well as within the 
7r~7r + 7r~ system. The isobar model decomposes the decay 
X~ — > 7r~7r + 7r~ into a chain of successive two-body decays: 
First the X~ with quantum numbers J PC decays into a di- 
pion resonance, the so-called isobar, and a bachelor pion 



bachelor 



The isobar has spin S and a relative orbital angu- 
lar momentum L with respect to 7r bachelol .- L and S couple 
to the X~ spin J. Second the isobar decays into tt + 7t~. This 
is illustrated in Fig. 6. In general the n~n + n~ system has 
isospin / > 1. Since there are no known flavor-exotic light- 
quark mesons and Pomeron exchange is assumed, we set 
/ = 1 . A state with odd number of pions has negative G- 
parity so that the value of the charge conjugation C is fixed 
by G = C(— l) 7 to be positive. 

The spin-density matrix pf ., is parameterized in terms 
of production amplitudes Vf : 



(5) 



N r 

r'l/ / j tr I'r 



where i = y 7>c M e [decay]L and e denote a particular partial- 
wave amplitude that is characterized by the isobar state, 
the decay orbital angular momentum L, the spin quantum 
numbers J PC of the X~ state, and by the spin projection 
M e of J. The amplitudes are constructed in the reflectivity 
basis [26], where M > and where the reflectivity e = +1 
describes the symmetry under reflection through the pro- 
duction plane. This basis is convenient, because, due to P 
parity conservation, amplitudes with different reflectivities 
do not interfere. Furthermore at high center-of-mass ener- 
gies -y^ m e reflectivity corresponds to the naturality of the 
exchanged Reggeon. The rank N r of the spin density ma- 
trix is set to two to account for the helicity flip and non-flip 
amplitudes at the baryon vertex assuming that the target 
nucleon stays intact. 
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The observed intensity is parameterized as a coherent 
and incoherent sum over the partial- wave amplitudes [26]: 



(6) I(r;m x ) = T 1 (T)Y J Yj X y ^^ (r;mz) 



Here t represents the five phase space coordinates that com- 
pletely describe the three-body kinematics. They are mea- 
sured for each event and are used to calculate the decay 
amplitudes i/^ which are constructed using non-relativistic 
Zemach tensors [27,28] and do not contain any free pa- 
rameters. The decay of the X~ is described in its rest sys- 
tem using the Gottfried- Jackson frame [29] with the z-axis 
along the beam particle direction and the {/-axis perpendic- 
ular to the production plane spanned by the beam and the 
recoil particle. The angular distribution of the isobar de- 
cay is defined in the canonical system obtained by pure 
Lorentz boost to the isobars rest frame. The production 
amplitudes Vf are complex numbers that are determined 
from extended maximum likelihood fits to the data per- 
formed in 40 MeV/c 2 wide bins in the three-pion invariant 
mass mx using a program that was originally developed at 
Illinois [30] and later developed at Protvino and Munich. 
This so-called "mass-independent" fit takes into account 
the overall acceptance r\ as a function of the phase space 
variables. It does not include any parameterization of the 
produced resonances X~ and assumes constant production 
strengths of the waves within the mx bin. The f depen- 
dence of the wave intensities was implemented by multi- 
plying the amplitudes with different f'-dependent functions 
of the form ff(t') oc e~ h< for M = and ff(t') oc t'e~ bt ' for 
M = 1, where the slope parameter b was determined from 
the data by fitting f' slices. 

The acceptance was estimated using a Monte Carlo 
simulation. Events were generated uniformly according to 
three-body phase space with beam kinematics and f' distri- 
bution taken from the real data and then processed through 
the detailed detector simulation and the full reconstruction 
and event selection chain. Figure 7 shows that COMPASS 
has an excellent acceptance of more than 50 % for n~ji + n~ 
final states even in the region mx > 2 GeV/c 2 . Moreover 
the acceptance shows only weak dependence on the po- 
lar angle of the isobar in the Gottfried-Jackson frame (cf. 
Fig. 8). This is important for the PWA, since in [20] it was 
shown that a drop of the acceptance towards cos 8qj = ±1 
in combinations with detector resolution effects may lead 
to significant leakage of intensity between waves. 

The PWA model includes five isobars [12]: The 
amplitudes of the p(770), / 2 (1270), and p 3 (1690) where 
described using relativistic Breit-Wigner line shape func- 
tions including Blatt-Weisskopf barrier penetration fac- 
tors [31]. The n + n~ S -wave is in general dominated by a 
broad /o(600) meson, a narrow /o(980), and several reso- 
nances above 1 GeV/c 2 . We use the parameterization 
from [33] which is based on the "M" solution from Au, 
Morgan, and Pennington [34], but with the /o(980) sub- 
tracted from the elastic nn amplitude and added as a sepa- 
rate Breit-Wigner resonance. 

Due to the angular momentum carried by the Pomeron 
that is exchanged in diffractive reactions, the incoming ii~ 
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(J p = 0~) can be excited to a state X~ with different J p , 
which is limited only by the conservation laws of strong 
interaction. In the construction of the PWA model waves 
with J <A and M < 1 were considered. The final wave set 
of the model consists of 41 partial waves plus one incoher- 
ent isotropic background wave, the so-called "flat" wave. 
Table 1 lists all waves. Mostly positive reflectivity waves, 
corresponding to natural parity exchange production, are 
needed to describe the data. The wave set in particular in- 
cludes the 2- + + [pn]F, 2- + \ + \pn]F, and 2- + l + [pn]P waves 
that in [20] were found to be crucial for describing the data 
and in addition made the exotic 7Ti(1600) signal disappear. 
Some of the waves have mx thresholds below which their 
intensity was fixed to zero. 

To ensure that the mass-independent fit found the global 
maximum of the likelihood, up to 100 independent fits were 
performed for each mass bin using random start parame- 
ters. Out of these fits the one with the largest likelihood 
was selected. In case two or more fits have ambiguous so- 
lutions that differ by less than one unit of likelihood, the 
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error of this bin was increased by the maximum difference 
of the solutions, which in the plots is indicated by thick 
green bars. 

After performing the mass-independent fit the mass de- 
pendence of the production amplitudes of a subset of six 
waves (printed in boldface in Table 1) was fit in the range 
from 0.80 to 2.32 GeV/c 2 to a model parameterized in 



Table 1. Wave set used for the mass-independent fit. The set con- 
sists of mostly positive reflectivity waves that correspond to nat- 
ural parity exchange. The six waves printed in boldface are used 
in the mass-dependent fit. 
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terms of Breit-Wigner amplitudes: 
(7) 

N r [resonances ] [resonances 

Plpn x ) = J] J] 4 r B1W J] c) lr BW ; (m x ) 

r V k J L I 



Here BWC^x) is the relativistic Breit-Wigner amplitude 
for resonance i: 



(8) BW i (m x ;M Q ,r ) = 



where r tot (m x ) = 2n c "' y r n (m x ) is the total mass-dependent 
width of the resonance including phase space factors and 
Blatt-Weisskopf barrier penetration factors B Ln (q n ) [31] for 
all A^decay decay channels: 



1 



m 2 -M? + ir tot (m x )M 



(9) 



r to t(m x ) = ^ r Q „ 



Mo q n B 2 Ln (q„) 
m x q „ B 2 Ln (q Qn ) 



with To = r tot (Mo). Here q„ represents the breakup mo- 
mentum of the particular two-body decay and qo„ = q n (Mo). 
If required by the data, for some waves an additional co- 
herent exponential background term of the form e~ aq was 
added. The waves selected for the mass-dependent fit ex- 
hibit either significant amplitudes or rapid phase motions 
in the 1.7 GeV/c 2 mass region. 



3.3 Results 

Figure 9 shows the most dominant wave, l ++ Q + [pn]S . The 
data points with statistical error bars are the result of the 
mass-independent fit. The wave intensity exhibits a broad 
structure around 1.2 GeV/c 2 which is the ai(1260). The 
continuous line represents the result of the mass-dependent 
^ 2 -fit of the intensities and phase differences of six waves 
(see Table 1). In this fit the peak is well described using 
a Breit-Wigner parameterization according to Bowler [35] 
and a small exponential background which could be caused 
by the Deck-effect [36]. 

The second most intense wave is the 2 ++ l + \pn]D wave 
which is depicted in Fig. 10. It is also the strongest M = 1 
wave and can be large only at high t' . The sharp peak of the 
a2(1320) is described by a Breit-Wigner, where the mass- 
dependent total width r tot (m x ) takes into account the two 
dominant decay modes of the 02(1320) to [pn]D and [qn]D 
with strength according to PDG [32]. In order to reproduce 
the high-mass tail of the wave and its interference with the 
l ++ + [pn]S and the 2~ + + [f2n]S waves, a second Breit- 
Wigner for the a 2 (1700) with its parameters fixed to the 
PDG values [32] was added to the fit. 

The third significant wave is the 2~ + + [/27r]S wave 
shown in Fig. 11. A Breit-Wigner corresponding to the 
7T2(1670) fits the data well. Its total width was parameter- 
ized assuming 60 % [f2n]S and 40 % [pn]P decays, ne- 
glecting contributions below 10 %. 

In order to study the sensitivity of the PWA to small 
signals, two low-intensity waves with known resonances 
were included in the mass-dependent fit: the A ++ \ + [pn]G 
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wave with the a 4 (2040) (see Fig. 12) and the (T + + [/ (980) 
tt]S wave with the tt(1800) (see Fig. 13). The A ++ \ + \pji]G 
wave peaks around 1.9 GeV/c 2 and is fit by a Breit-Wigner 
amplitude using a constant width in its denominator, be- 
cause no branching fractions are known for the 04(2040). 
The resulting 04 mass of 1 884 ± 13(stat.)+f (syst.) is lower 
than the PDG average [32], but has a large systematic er- 
ror. Even fainter in intensity is the tt(1800) peak in the 0~ + 
+ [/o(980)7r]S wave, but it exhibits a clear phase motion 
with respect to the 7^(1670) in the 2~ + + [/27r]S wave (not 
shown here). A Breit-Wigner with constant total width on 
top of a background describes the data. 

Figure 14 shows the intensity of the spin-exotic 1" + 
l + [pn]P wave. The distribution exhibits a broad bump at 
1.7 GeV/c 2 and a shoulder at lower masses. The phase dif- 
ference between this wave and the l ++ + [pn]S wave, as 
depicted in Fig. 15, exhibits a rising phase motion around 
1.7 GeV/c 2 , in the tail region of the ai(1260), thus indi- 
cating a resonant behavior of the spin-exotic wave. This 
is supported by the flat phase difference with respect to 
the 2~ + + [/27r]S wave as shown in Fig. 16 which can be 
explained by a 1~ + resonance with mass and width simi- 
lar to that of the 7^(1760) which dominates the 2~ + wave. 
The intensity and the phase differences of the 1~ + ampli- 
tude are fit by a Breit-Wigner with constant total width 
on top of an exponential background that describes the 
low-mass shoulder. The extracted resonance parameters of 
M = 1 660 ± 10(stat.)^ 4 (syst.) MeV/c 2 and T = 269 ± 
21(stat.)^(syst.) MeV/c 2 are in agreement with the dis- 
puted 7Tl(1600). 

Table 2 summarizes the resonance parameters extracted 
by the mass-dependent fit which are in good agreement 
with the PDG values [32]. The table also gives the intensi- 
ties of the resonant part of the particular waves integrated 
over the mass range from 0.80 to 2.32 GeV/c 2 and nor- 
malized to the total intensity of the waves in the mass- 
dependent fit, which corresponds to 38.7(2) % of the total 
acceptance-corrected data sample in the mass range. 



xlO 3 




Mass of 7i"7i"7i + System (GeV/c 2 ) 

Fig. 9. Intensity of the l ++ + \pn]S partial wave. The wave ex- 
hibits a broad peak from the ai(1260). The continuous line shows 
the result of the mass-dependent fit. 
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Fig. 10. Intensity of the 2 ++ l + [pn]D partial wave. The wave is 
dominated by a narrow a2(1320) peak. The continuous line shows 
the result of the mass-dependent fit. 




Fig. 11. Intensity of the 2~ + + [f 2 n]S partial wave. The wave ex- 
hibits a broad 7^(1670) peak. The continuous line represents the 
result of the mass-dependent fit. 



Several studies were performed to test the stability of 
the 1~ + wave with respect to various assumptions made in 
the analysis, e.g. by adding and removing waves, varying 
cuts or initial values of the fit parameters, and shifting the 
mass binning in the mass-independent fit. In addition the 
dependence of the fit on the rank N r was studied. The fit 
with N r — 1 shows an increase of the flat background wave 
from 5.8 % for N r = 2 to 19 % with respect to the total 
acceptance-corrected data sample in the mass range from 
0.5 to 2.5 GeV/c 2 . For N r = 3 this value drops to 1.2%. 
This fit, however, exhibits larger bin-to-bin fluctuations. In 
both cases the 1~ + intensity was not significantly altered in 
the region between 1.5 and 1.8 GeV/c 2 . 

In [20] it was shown that an inhomogeneous accep- 
tance not completely taken into account by the Monte Carlo 
simulation or an incomplete wave set may lead to leakage 
of intensity of non-exotic waves into the 1~ + 1 + wave. This 
effect was studied by generating Monte Carlo events using 
the parameters of the 16 most dominant waves excluding 
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Table 2. Preliminary resonance parameters and intensities of the specified decay channels for the six waves included in the mass- 
dependent fit. The first uncertainty corresponds to the statistical error, the second one to the systematic error. 



State 


Mass 
[MeV/c 2 ] 


Width 

[MeV/c 2 ] 


Intensity 

[%] 


Channel 
/ pc Af e [decay]L 


a, (1260) 


1255 ± 


6% 


367 ± 9 + _f 5 


67 ± 3 + _ 4 20 


l ++ + [pn]S 


a 2 (1320) 


1321 ± 




110 ± 2 + _ 2 15 


19.2 ± 0.6+oj 


2 ++ l + [pn]D 


7r 2 (1670) 


1 658 ± 


a+24 
J -8 


271 ± 9 + _H 


10.0 ± 0.4^ 


2- + 0\f 2 x]S 


a 4 (2040) 


1 885 ± 


13+f 


294 ± 25^9 


i.o ± 0.3+g ; 


4 ++ l + [/OT]G 


tt(1800) 


1 785 ± 


q+12 


208 ± 22+2' 


0.8 ±0.1+°] 


Q- + + [f n]F 


/ri(1600) 


1660± 


10& 


269±21 + g 


1.7 ±0.2+° ? 
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Fig. 12. Intensity of the 4 ++ 1 + [/ot]G partial wave. The broad peak 
around 1.9 GeV/c 2 is attributed to the a 4 (2040). The continuous 
line shows the result of the mass-dependent fit. 
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ture 
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13. Intensity of the 0" + + [/ (980)ff]5 partial wave. The struc- 
around 1.8 GeV/c 2 is attributed to the 7r(1800). The continu- 
line shows the result of the mass-dependent fit. 



the 1 + wave. The Monte Carlo data were then processed 
by the same PWA used for the real data. The amount of 
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Fig. 14. Intensity of the spin-exotic l~ + \+[pn]P partial wave. The 
distribution exhibits a broad bump around 1.7 GeV/c 2 . The red 
curve shows the result of the mass-dependent fit with one Breit- 
Wigner for the n t ( 1 600) (blue curve) on top of a background (pur- 
ple curve). 



intensity leaking from the 16 major waves into the 1~ + 1 + 
[pn]P wave was found to be below 5 % of the 1~ + inten- 
sity and is therefore negligible. This is mainly due to the 
flat angular acceptance and the better resolution of COM- 
PASS, both well reproduced by the simulations. 

In the mass-dependent fit the low mass shoulder in the 
1~ + 1 + |/ot]P wave is accounted for by a non-resonant expo- 
nential background, which is possibly caused by a Deck- 
like effect [36]. In an alternative attempt to describe the 
shoulder an additional Breit-Wigner for the 7^ (1400) with 
its parameters fixed to the PDG values [32] was included 
in the mass-dependent fit. This did not significantly alter 
the intensity or the phase difference of any of the waves 
in the mass-dependent fit, but slightly shifted the 7Ti(1600) 
towards smaller masses which was included in the system- 
atic error. 

Also using different parameterizations for the nn S- 
wave and the p(770) and replacing the Zemach tensors by 
rotation functions with relativistic corrections [37] did not 
significantly change the result. This is also true for tests, 
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Fig. 15. Phase difference of the \~ + \ + \pnW and the l ++ + |p^]5 
partial waves. A clear phase motion is seen in the region around 
1.7 GeV/c 2 which can be explained by a 7T[(1600) resonance in- 
terfering with the tail of the ai(1260). The continuous line shows 
the result of the mass-dependent fit. 



150 



$ 100 



Oh 



50 

-50 
-100 
-150 
-200 



Jt Pb — > 7T"7t"Ji; + Pb 

- O.Kf < 1.0 GeV 2 /c 2 



COMPASS 2004 



ii 



A^)(l" + l + pjtP-2" + + f 2 7lS) 

t 1 1 I 1 1 I L I [ 1 I I I 1 1 I I 1 1 t I 



0.6 0.8 1 1.2 1.4 1.6 1. 



2 2.2 2.4 



Mass of n'ii'n + System (GeV/cr) 

Fig. 16. Phase difference of the r + l + [p^]P and the 2- + + [f 2 n]S 
partial waves. The phase difference is virtually flat which can 
be explained by a 7Ti(1600) resonance with resonance parame- 
ters similar to the 7r 2 (1670). The continuous line shows the result 
of the mass-dependent fit. 



where the PWA was performed without the t' -dependent 
factors in the partial-wave amplitudes. 

It cannot be ruled out, however, that the 7Ti(1600) sig- 
nal is an artifact of the simplifying assumptions made in 
the PWA model. Final state interactions in particular the 
Deck-mechanism were found to be relevant in the 
system [38,39] and might cause spurious signals. 



mostly on liquid hydrogen targets. These data will make it 
possible to study a number of channels with unprecedented 
precision. 

Several analyses have started: Figure 17 shows the 
7r~7r + 7r~ invariant mass spectrum after event selection. The 
tremendous boost in statistics will enable COMPASS to 
extract precise resonance parameters, to study the f' -depen- 
dence of resonance production, and to hopefully settle the 
7Ti(1600) issue in the diffractive ji~n + n~ channel. This anal- 
ysis will be accompanied by a PWA of the n~7i {) n a final 
state which will provide important internal cross checks of 
our results. Furthermore the PWA will be extended to the 
region of t' < 0.1 (GeV/c) 2 . Data sets with different targets 
will give access to possible nuclear effects in diffractive 
production. 
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Fig. 17. Invariant mass distribution n n + n final states diffrac- 
tively produced by a n~ beam on a liquid hydrogen target; based 
on 21 % of the 2008 data set. 



In addition COMPASS will search for scalar and tensor 
glueballs in centrally produced n + n~ji + n~ systems. Fig- 
ure 18 depicts the 7r + 7r~7r + 7r~ invariant mass with and with- 
out selecting centrally produced event via the cut x F (n Ust ) > 
0.7. A clear signal from the /i (1285) is seen. 

The f\7i and maybe even the b\n channel, into which 
the 1~ + hybrid is predicted to preferentially decay, can also 
be studied in charged five pion final states produced in 
ttA — > 7T~7r + 7T~7T + n~A. Furthermore this channel gives ac- 
cess to higher masses above 2 GeV/c 2 , where the PDG lists 
many states that need confirmation [32]. 

Finally COMPASS will search for possible hybrids and 
glueballs in kaonic final states produced in reactions like 
TCp -> n K + K p, n-K Q s K Q s p, or tT£<!#°/j. Moreover the 
CEDAR detectors enable COMPASS to measure diffrac- 
tion of kaon beam, e.g. in K p — > K n + n'p. 



3.4 Run 2008/9 Data 

The analysis of the data from the 2004 pilot run marked 
only the beginning of the hadron spectroscopy program at 
COMPASS. During the years 2008 and 2009 COMPASS 
collected very large diffractive and central production data 
sets using 190 GeV/c positive and negative hadron beams 



4 Summary 

COMPASS has started its hadron spectroscopy program 
with the main goal to precisely explore the light-meson 
sector in order to search for states with gluonic degrees of 
freedom like hybrids and glueballs and to settle the proper- 
ties of controversial states like the 7Ti(1400), 7Ti(1600), and 
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Fig. 18. Invariant mass distribution of the n + n~n + n~ system pro- 
duced in n~p — > 7Tf ast (;r + 7r~7r + 7r~)p recoi | with and without the cut 
Xjr(^ fesl ) > 0.7; based on 13 % of the 2008 data set. 

the /o(1500). At COMPASS two production mechanisms, 
diffractive dissociation and central production, can be mea- 
sured in parallel using n, K, or p beams. The spectrometer 
provides high acceptance and excellent resolution for final 
states containing both charged and neutral particles. 

In a short first pilot run in 2004 COMPASS acquired a 
diffractive data set using a 190 GeV/c 7T~ beam scattered off 
a lead target. A partial wave analysis of the n~n^ii~ final 
state in the region of t' e [0.1, 1.0] (GeV/c) 2 shows sev- 
eral of the well-known resonances expected in this chan- 
nel. In addition we observe a spin-exotic J PC = 1~ + state 
with significant intensity in the pn decay channel in nat- 
ural parity exchange at a mass of 1.66 GeV/c 2 . The mass 
dependence of its phase and its resonance parameters are 
consistent with the disputed ;ri(1600). 

Based on the experience of the pilot run the spectrom- 
eter was upgraded and improved. In particular a detector 
to measure the target recoil was added to the setup. Dur- 
ing dedicated runs in 2008 and 2009 COMPASS accumu- 
lated diffractive and central production data sets of un- 
precedented statistics using 190 GeV/c hadron beams. The 
bulk of the data were taken with a liquid hydrogen target, 
but also nickel, tungsten, and lead targets were measured. 
First analyses of diffractive and central production reac- 
tions have started. 
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